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We study how local fluctuations in the initial states ol relativistic heavy-ion collisions manifest 
themselves in the correlations between different orders of harmonic moments of the density proflles, 
particularly those involving only odd harmonics which purely arise from initial state fluctuations. 
We flnd the strengths of those correlations are sensitive to the number of hot and cold spots in the 
initial states. Hydrodynamic evolution of the flreball translates initial state geometric anisotropies as 
well as their correlations into final state momentum anisotropies and correlations. We conclude that 
the measurement of the correlations between different harmonic moments of final state azimuthal 
distribution can be employed to quantify the inhomogeneity of the initial density profiles such as 
the population of hot and cold spots that are produced in high energy nuclear collisions. 



The creation of hot, dense quark-gluon plasma (QGP) 
can be achieved by colliding two heavy nuclei at ultra- 
relativistic energies, such as those at the Relativistic 
Heavy-Ion Collider (RHIC) and the Large Hadron Col- 
lider (LHC). The dynamical evolution of the QGP pro- 
duced in these energetic collisions, particularly the large 
collective flow observed at RHIC and the LHC [l|4j , can 
be well simulated by relativistic hydrodynamics [1]. Be- 
ing the hydrodynamic response to the pressure gradients, 
the collective flow exhibited by such highly exited QCD 
matter is azimuthally anisotropic in the plane transverse 
to the beam axis in the collisions with nonzero impact 
parameter. The anisotropy of the flow is usually quanti- 
fied by the Fourier expansion coefficients 7j„ of the final 
state momentum distribution in the transverse plane Q . 
There has been extensive study of elliptic fiow V2 as a 
function of various quantities, mostly aiming for a quan- 
titative extraction of the transport properties, e.g., the 
shear viscosity to entropy ratio, of the produced hot QCD 
matter 

Early comparisons to the measured elliptic flow V2 us- 
ing relativistic hydrodynamic simulation employed a set 
of smooth, event-averaged initial conditions which are 
specified by the event-averaged geometry of the collision 
zone, usually in terms of energy or entropy density at 
initial times. However, the outcome from each nucleus- 
nucleus collision fiuctuates event- by- event due to quan- 
tum fluctuations, such as the positions of nucleons or 
color charges inside the colliding nuclei [lol - [T3| . One 
consequence of such fluctuations in the initial states is 
the finite eccentricity and elliptic flow even in the colli- 
sions with almost zero impact parameter [l^ . Another 
significant feature of the initial state fiuctuations is the 
presence of odd harmonic moments in the initial geo- 
metric anisotropy and final momentum anisotropy |15l - 
[28j . in addition to even harmonic moments which are 
just the refiection of the symmetry of the overlap re- 
gion between two colliding nuclei at finite impact pa- 



rameters. The fiuctuating non-smooth initial conditions 
have been invoked to explanin the double-hump struc- 
ture and near-side ridge phenomenon in two-particle cor- 
relation measurements [l^ IsoMsil . The measurements of 
azimuthal anisotropy of the final momentum distribution 
from RHIC and the LHC have shown a prominent third 
harmonic flow and other odd harmonic flows, in ad- 
dition to even harmonic flow coefficients (35l-l37| . Signifi- 
cant attention is now devoted to analyzing the dynamical 
evolution of initial state fiuctuations and their hydrody- 
namic responses and translation to final harmonic fiows 
B [13, HO, m, HI- The goal is to achieve a quantitative 
understanding of the expansion dynamics of the fireball 
produced in these energetic nucleus-nucleus collisions. 

In this work, we study another important aspect of 
initial state fluctuations: the inhomogeneity of the ini- 
tial density profile. The local fiuctuations (hot or cold 
spots) not only generate geometric anisotropy in terms 
of harmonic moments and their fluctuations, different or- 
ders of harmonic moments are actually correlated with 
each other. Some of this type of correlations have been 
investigated in earlier studies [l9l . [23l [27| . Our objective 
is to study how the local fluctuations on top of event- 
avcragcd profiles manifest in the correlations between 
different harmonic moments of the initial state geomet- 
ric anisotropy. Due to hydrodynamic evolution of the 
fireball, these correlations will be translated to the cor- 
relations between final state particles. Our particular 
interests are focused on those correlations involving only 
odd harmonic moments which purely stem from the lo- 
cal fiuctuations in the initial states. We show that the 
correlations between different harmonic moments can be 
employed to study the degree of inhomogeneity of the 
initial density profile. Especially, the magnitudes of the 
correlations strongly depend on the number of hot / cold 
spots present in the initial states. Once measured by ex- 
periments, these correlations may infer much information 
about the initial state fluctuations, such as the popula- 
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tion of hot/cold spots inside the hot, dense QGP pro- 
duced in ultra-relativistic nuclear collisions. 

For a given initial energy or entropy density profile 
/(r) = /(r, (/)) in the transverse plane, we may quantify 
the anisotropy of the profile by defining eccentricities e„. 



{r'"cos(n0-n$„)}/{^"} 
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Here we use curly brackets {• • • } to represent the average 
of density profile within a given event and angle brackets 
(• • • ) for the average of quantities over many events (the 
ensemble average). The initial spatial event plane angle 
<i>„ for the n-th harmonic moment is define as 

$„ = - arctan [{r™ sin(?i(/))}/{r" cos(n(/))}] (2) 
n 

where (j) =^ arctan(?//x) represent the polar angle for the 
point {x,y). Here we take the exponent in r™ to be 
m = n for n > 2 [20j and to = 3 for the dipole asymme- 
try n = 1 following Teaney and Yan [23| , who introduced 
a cumulant expansion to parameterize the initial condi- 
tions. Gardim et al. [S^ have argued that the choice 
adopted here serves a better estimator than m = 2 for 
the flow harmonics u„ with n > 2. 

Being the hydrodynamic response to the initial state 
spatial anisotropy e„, the momentum anisotropy param- 
eters Vn (flow harmonics) of the final state particle az- 
imuthal angle {ij}) distribution are defined as 



Vn = {cos(ri-0 - n*„)} 



(3) 



where 4'„ is the final (momentum) event plane angle. In 
an event-by-event hydrodynamic analysis with fiuctuat- 
ing initial conditions, the final momentum event plane 
is strongly correlated to the initial spatial event plane $„ 
[TtI [25l.l39j. Assuming linear hydrodynamic responses of 
anisotropic flows Vn to the initial spatial anisotropy e„, 
one may set '^n = $n + tt/u. Note there may exist non- 
linear hydrodynamic response of w„ to e„, e.g., V4 can be 
developed from £4 as well as £2- 

In relativistic nuclear collisions, the initial states fluc- 
tuate from one event to another due to local fluctuations. 
Thus for a given event we may write down the initial den- 
sity profile f{r) in the transverse plane as the sum of an 
event-averaged profile foir) and an event-by-event fluc- 
tuating one Sf{r), 



/(r1-/o(r) + <5/(rO 
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By definition, (/(f)) = /o(f) and {6f{f)) 0. The event- 
averaged profile /o('^ is usually specified in a particular 
model; here we use the optical Glauber model [i^ |4l| , 
with the initial density being proportional to the combi- 
nation of binary collision density and participant nucleon 
density [42|: /o(f) oc [ancoii(f) + (1 - a)np^rt{r)/2]. For 
the simulation of the collisions of two lead-lead collisions 
at ^/sNN = 2.76 TeV at the LHC, we use Woods-Saxon 
profile for the the nuclear density, with the radius pa- 
rameter taken to be /J^g = 6.62 fm and the diffusion 
parameter d = 0.546 fm. The inelastic nucleon-nucleon 



cross section is taken to be crjvAr = 64 mb. The param- 
eter a is chosen as 0.11 for a good description of the 
centrality-dependence of the charged hadron multiplicity 
distribution 4^. 

The harmonic moments e„ and $„ can be expressed 
using event-averaged and fluctuating profiles. They are 
particularly simple for odd harmonic moments: 



— arctan [{?-" 
n 



' sin(n0)}5/{r'" cos(n0)}5] 
en = {r" cosincj) - n<^>n)}s/{{r"'}o + {r™}^) (5) 

where the subscripts in {• • • }o and {• • ■ represent the 
average over the profiles /o(f) and 6f(r), respectively. 
We see that for odd harmonic moments, the event-plane 
angles arc completely determined by the fluctuating 
part of the profile Sf{r), and the event-averaged profile 
fo^r) only contributes to the overall normalization for the 
eccentricities e„. 

For given event-by-event initial conditions including 
various fluctuations, we may obtain the event distribu- 
tion of the fluctuating part of the profile 6f{r). In this 
work, we model the fluctuating part 6f{r) for a given 
event by a combination of N local fluctuations. 
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Here, fi = {ri,(j)i) is the location of the fluctuation with 
radius of Rh (taken to be the same for all fluctuations). 
The sign factor Si — ±1 represents whether the local 
fluctuation is positive or negative (hot or cold spot) and 
Aii{si) represents the amplitude of the fluctuation. Both 
hot and cold spots are generally called hot spots in what 
follows. The local fluctuations are assumed to be small 
perturbations on top of the event-averaged profile, i.e., 
Ah{±l)N < J (frfo{f). The probabilities for hot and 
cold spots are related to their amplitudes by Pg. (s^ = 
l)Ahil) = PsAs^ = -l)Ahi-l) to ensure (5/(f)) = 0. 

The description of initial state fluctuations with hot 
and cold spots is the approximation of the realistic event 
distribution P{Sf) by two spikes located at A/i(l)/(7ri?^) 
and — A/i(— l)/(7ri?^) with the heights being Psi(l) and 
Ps;(— 1), in addition to a spike at zero (whose height, 
denoted as Psi(O), is not relevant to the following discus- 
sion). The heights of two spikes are usually not equal 
since the event distribution P{Sf) in realistic initial con- 
ditions is not symmetric around zero. 

As mentioned earlier, the presence of hot spots gener- 
ate not only the spatial anisotropy in terms of harmonic 
moments e„, but the correlations between different or- 
ders of harmonic moments as the hot spots arc coherent 
combinations of all orders of harmonic moments. In par- 
ticular, the event plane angles $„ of different orders of 
harmonics are correlated to each other. Such correlations 
can be quantified by the functions C^''-*(ni, n2 - ■ - rik), 

C'^'') = (cos[rii$„, -I- n2$„, + • • • + nfc<I>„J) (7) 

Here we are interested in the correlation functions with 
the combinations of indices being ni + n2 + ■ ■ ■ + Uk = 
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as these are not dependent on the particular choices of 
the coordinate. For relativistic heavy-ion colhsions, the 
above correlation functions arc completely determined by 
the early time dynamics of the collisions. Hydrodynamic 
evolution of the fireball will translate these correlations 
into the correlations in the final states. With the substi- 
tution of = $„ -|- 7r/n, wc may cast the above corre- 
lation functions in terms of final momentum event plane 
angles, 

CW = (_l)fe(cos[ni*„, + n2*„, + • • • + nfe*„,]) (8) 

This simplified substitution will be modified by the non- 
linearity of the hydrodynamic equations, however, we 
expect that any such correlations in the initial state 
will survive the hydrodynamic evolution and manifest 
themselves in the final state, especially when the fluc- 
tuations have small amplitude. Of particular interests 
are these involving only odd harmonic moments as they 
are purely from initial state fluctuations and less af- 
fected by the collision geometry. For illustration pur- 
pose, here we consider these combinations involving only 
the first three odd harmonic moments, C^'^''{1, f, 1,-3), 
C(4)(l,l,3,-5) and C^^) (3, 3, -1, -5). 

To investigate how the hot spots on top of event- 
averaged profile affect the correlations between differ- 
ent harmonic moments, we need to model the spatial 
distribution of hot spots in the initial states. Here we 
take two typical distributions of hot spots for compari- 
son purpose. The first one is the hard sphere model in 
which hot spots are sampled according to a uniform dis- 
tribution in the overlap region of two hard sphere nuclei: 

Pffs(r) cx e{RHs - |r - 6/2|)0(i?Hs - \r + b/2\). The 
other is taken as the two-component Glauber model in 
which hot spots are more distributed towards the center 
compared to the hard sphere model. 

In Fig.[Tl we show the event-averaged correlation func- 
tions as a function of the number of hot spots in the 
initial states. The panels (a, b) show the results for the 
hard-sphere model of hot spot distribution, and (c, d) for 
the Glauber model. The panels (a, c) show the results 
for the isotropic event distribution of hot spots in the 
transverse plane (6 = 0), while (b, d) for anisotropic dis- 
tribution (6 7^ 0) . Here for anisotropic hot spot distribu- 
tions, we choose the impact parameter 6 = 12 fm between 
two colliding nuclei (hard spheres or Woods- Saxon pro- 
files); the hard sphere radius Rhs is tuned to give similar 
values of eccentricity €2 as the Glauber model when cal- 
culated with the event distributions of hot spots. The 
radius Rh of the hot spots is taken to be 0.25 fm. The 
increase of hot spot radius tends to decrease the cor- 
relations between different harmonic moments since hot 
spots are more evenly distributed in a given event. Such 
effect is small for large system size or small number of hot 
spots. The relative probabilities of hot and cold spots are 
taken as Psi(l) = 35% and Ps.{—1) = 65% by compar- 
ing to the values of the correlation functions C^^\ e.g., 
C(^'(2,3,-5) « 0.1 from Monte-Carlo Glauber modeling 
of lead-lead collisions at 6 = fm and our model with 



about a hundred hot spots. 
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FIG. 1: (Color online) The correlation functions C''*' as a 
function of the number of hot spots in the initial states. 
The panels (a, b) are for the hard sphere and (c, d) for the 
Glauber model; the panels (a, c) are for isotropic and (b, d) 
for anisotropic event distributions of hot spots. 

The effects of the hot spots on these correlation func- 
tions are clearly seen from these plots. If there is only one 
hot spot in the initial state, all correlation functions are 
equal to unity since all event plane angles are aligned or 
anti-aligned to each other, the strongest correlation that 
can be achieved. With increasing number of hot spots, 
the interference between different hot spots decreases the 
correlations between different harmonic moments. The 
geometric effect on the correlation functions can be seen 
by comparing (a) and (b) [or (c) and (d)] panels. When 
the initial state contains few hot spots in, such effect is 
small since the hot spots are barely correlated with the 
geometry in each event. But when the number of hot 
spots is large, all the event plane angles are strongly bi- 
ased by the geometry in each event, thereby enhancing 
the correlations between different harmonic moments. 

Due to different spatial (r) distribution of hot spots, 
there exists some difference in the correlation functions 
between these two models. In central collisions, the cor- 
relations are slightly stronger in the Glauber model dis- 
tribution of hot spots. This is due to the fact that the 
hot spots are more distributed toward the center, thus 
fewer hot spots lie close to the edge of the overlap re- 
gion. This also reduces the dependence on the geometry 
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of the hot spot distribution, i.e., smaller enhancement of 
the correlations from isotropic distribution of hot spots 
to anisotropic one in the Glauber model compared to the 
hard sphere model. 
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FIG. 2: (Color online) The correlation functions C^^' from 
Monte-Carlo Glauber modeling of the initial conditions for 
Pb+Pb collisions at y/sJ^N = 2.76 TeV at the LHC. 

Next we calculate the correlations between different 
harmonic moments for realistic simulations of lead-lead 
collisions at the LHC energy of ^^snn = 2.76 TeV us- 
ing the Montc-Carlo Glauber model following Rcf. [20| . 
where the fluctuations of nucleon positions inside two 
colliding nuclei as well as the fluctuations arising from 
individual nucleon-nucleon collisions are included. The 
results for the correlation functions C^^^ are plotted as 



a function of impact parameter b in Fig. [2] Here we 
take the initial condition with free-streaming of 0.25 fm 
to account for the evolution of hot spots during the pre- 
equilibrium stage. One can see that the correlations be- 
tween different harmonic moments increase from central 
collisions to peripheral collisions, indicating a decrease in 
the number of hot spots. Compared with the results from 
Fig. [TJ we obtain just a few hot spots in very peripheral 
collisions and up to about one hundred hot spots in most 
central collisions. We note that the number of hot spots 
may be different in other initial condition models, such 
as the Color Glass Condensate model [43^46} , where the 
number of fluctuations is essentially determined by the 
transverse size of the color flux tubes and the overlap 
area of the colliding nuclei. 

In summary, we have studied the effect of local fluctu- 
ations in the initial states on the geometric anisotropy by 
considering the correlations between different harmonic 
moments. We find that the strength of such correlations 
strongly depends on the number of hot spots produced 
at initial times. Such correlations in the initial states can 
be directly related to the correlations in the final states 
owing to the hydrodynamic response to the pressure gra- 
dient of the fireball. Thus the measurement of the corre- 
lations between four different harmonic moments, if ex- 
perimentally feasible, can determine how many hot spots 
are inside the quark-gluon plasma created in relativistic 
heavy-ion collisions. 
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